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Finding a suitably growing length scale that increases in tandem with the immense viscous slowdown of 
supercooled liquids is an open problem associated with the glass transition. Here, we define and demonstrate 
the existence of one such length scale which may be experimentally verifiable. This is the length scale over 
which external shear perturbations appreciably penetrate into a liquid as the glass transition is approached. 
We provide simulation based evidence of its existence, and its growth by at least an order of magnitude, by 
using molecular dynamics simulations of NiZr 2 , a good fragile glass former. On the probed timescale, upon 
approaching the glass transition temperature,from above, this length scale, is also shown to be consistent 

with Ising-like scaling, f oc j , with y ^ 0.7. Furthermore, we demonstrate the possible scaling of f about 

the temperature at which super-Arrhenius growth of viscosity, and a marked growth of the penetration depth, 
sets in. Our simulation results suggest that upon supercooling, marked initial increase of the shear penetration 
depth in fiuids may occur in tandem with the breakdown of the Stokes-Einstein relation. 

PACS numbers: 75.10.Jm, 75.10.Kt, 75.40.-s, 75.40.Gb 


Introduction. When a liquid is cooled sufficiently quickly 
to temperatures well below its melting temperature, nucle- 
ation is avoided and the transition to the crystalline state, pos¬ 
sessing both extended long-range structural order and abso¬ 
lute minimum free energy, is bypassed. A liquid maintained 
beneath its melting temperature exists in metastable equilib¬ 
rium, and is said to be supercooled. A supercooled liquid 
lacks the long range structural order characterstic of the under¬ 
lying crystalline ground state, instead maintaining the amor¬ 
phous atomic arrangement typical of a liquid. As the tem¬ 
perature of the supercooled liquid is lowered further, the vis¬ 
cosity (and relaxation time) increases dramatically, by up to 
some 14 decades over a temperature range as small as 100 K. 
Eventually, a temperature, Tg, is reached at which the viscos¬ 
ity (and hence, relaxation time) is so large (>10^^ Poise/100 
s) that structual rearrangments cease to take place on any rea¬ 
sonable timescale, and the liquid behaves rigidly in response 
to fluctuations and perturbations. By definition, the liquid is 
then out of equilibrium, and this is deemed the glass transi¬ 
tion. The ’transition’ occuring at Tg, is in fact not a thermo¬ 
dynamic transition, but instead a kinetic crossover. There is 
no thermodynamic driving force (energy saving) associated 
with Tg, and a structural rearrangment and associated symme¬ 
try breaking is apparently absent. In addition to the smooth 
emergence of rigidity at Tg, the glass transition is accompa¬ 
nied by a rich phenomenology and wide range of interesting 
features that cannot be enumerated here, but are discussed in 
a variety of exceptional reviews, e.g. EM- 

The two most puzzling aspects of the glass transition are the 
onset of structural rigidity without apparent long range struc¬ 
tural order (with associated long-time, non-zero shear modu¬ 
lus), and the dramatic, faster than Arrhenius increase of the 


* weingartner.n.b@wustLedu 
zohar©wuphys.wustl.edu 


viscosity/relaxation time found in the so-called fragile glass 
formers HZlIEl. Both features seem to call for, and likely re¬ 
quire, the existence of a growing length scale, intimately con¬ 
nected to the propagation of some form of amorphous order 
or increasingly cohesive, extended network. In fact, while 
the notion of a growing activation energy barrier is clearly 
tied to cooperative motion, even more fundamentally, sim¬ 
ple intuitive reasoning suggests that a dramatically growing 
(or diverging) timescale to relaxation should be coupled to 
a similarly increasing (and possibly diverging) length scale. 
Recently, rigorous bounds mandating the existence of a con¬ 
comitant growth of spatial length scale with relaxation time 
have been proven to exist |0| . 

The notion of a growing length scale underlying the dy¬ 
namic slowdown of the glass transition is a principle feature 
of most theories of glass formation, such as that of Adam and 
Gibbs, Random First Order Transitions, Mode-Coupling The¬ 
ories, Kinetically Constrained Models, and others CHSlIIo]- 
[TSl. Some of these theories predict an underlying phase tran¬ 
sition at a temperature below Tg, with the glass transition serv¬ 
ing as a kinetic ’’ghost” preceding the actual thermodynamic 
change. Others posit that there is no true thermodynamic 
transition besides the melting/freezing transition, and that the 
length scale corresponds to a geometrically arrested structural 
ordering which is still capable of bringing about rigidity. As 
such, the quest to find physical, verifiable, and suitably in¬ 
creasing length scales has been underway for decades. Many 
proposals for appropriate length scales have been made in¬ 
cluding those associated with liquid-like defects, the lowest 
eigenvalues of the relevant Hessian matrix for a system, var¬ 
ious point-to-set lengths, elasticity lengths |[T9H34ll . and dy¬ 
namical heterogeneity lengths |[35ti43]| and computer vision 
methods to ascertain both static and dynamic length scales 
1441 . ’’Hybrid” correlation length scales, that have mixed 
static/dyanmical charactersitics, have also been found, as in 
1451 . Each of the previously proposed length scales is excep- 
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tionally interesting in their own right (and perhaps many can 
eventually be found arise from the same underlying mecha¬ 
nism), but many display the same drawbacks. Previous nu¬ 
merical and experimental work has shown that these length 
scales evade experimental verification, and/or do not display 
an exceptional growth upon approach to Tg. For instance, in 
the case of the lengthscale investigated in EH, the behav¬ 
ior of the lengthscale and underlying physics bears a passing 
resemblance to the investigation done in this work, but that 
lengthscale, which arises in response to internal perturbations 
associated with thermal fluctuations, requires knowledge of 
individual particle displacements, making it difficult to exper¬ 
imentally detect. It is natural to expect that propagating amor¬ 
phous order should be able to be revealed experimentally, and 
our proposed correlation length has the benefit of being read¬ 
ily measurable with methods beyond scattering experiments. 

It is also worth noting that in 1461 . the authors found ev¬ 
idence of a decreasing correlation length upon cooling to¬ 
ward the glass transition temperature, Tg, in kinetically strong 
glassforming liquids. This was found to be indicative of a 
lambda transition in the vicinity of the glass transition. It was 
further suggested that for kinetically fragile liquids, the be¬ 
havior would be the opposite, with the correlation length in¬ 
creasing upon approach to Tg. While this general behavior 
is consistent with our findings, we found no evidence sug¬ 
gestive of a lambda transition in this system, and further, the 
correlation length investigated in ll46l . was purely dynamical 
in nature. 

In light of the above discussion, and based upon sugges¬ 
tions made in previous theoretical work Eaim, we perform 
molecular dynamics simulations of NiZr 2 , an excellent repre¬ 
sentative of a fragile glass ||48]| . We provide evidence for the 
rapid increase of the shear penetration depth, defined as the 
length over which a supercooled liquid rigidly responds to ex¬ 
ternally imposed forces. We find that the near divergence of 
the penetration depth as the system becomes glassy is not far 
off the mark of Ising-like scaling. 

The Shear Penetration Depth. 

In ordinary critical phenomena, the correlation length scale 
is defined as the typical spatial extent of a fluctuation of the 
thermodynamically relevant order parameter. It can also be 
interpreted as the average length over which a perturbation 
by the appropriate conjugate ’’generalized force” will appre¬ 
ciably propagate. For example, in the Ising Model, the cor¬ 
relation length corresponds to fluctuations in the typical size 
of magnetic domains (the order pameter is the magnetization, 
M), and the also corresponds to the distance over which an 
applied magnetic field, ^ (the conjugate force) will influence 
the system. As it is known that in crstalline solids, the rigidity 
is due to long range order, we can apply this idea to the glass 
transition problem. In this case the ordering should be short 
range at temperatures just below melting, and grow as temper¬ 
ature as lowered. We can quantify this by subjecting the liquid 
to a shear perturbation on the boundary, and tracking how it 
penetrates the liquid transverse to the applied stress. Liquids, 
by definition, are capable of rearrangement to dissipate shear 
stress; a shear force applied to the top of a liquid will only 
propagate appreciably through a finite number of layers be- 





FIG. 1. (Color Online) Representation of the proposed response of 
general supercooled fluid systems. The solid lines represent the orig¬ 
inal box shape before perturbation. The dashed regions represent 
the successive layers that respond to the perturbation at temperatures 
above and around Tg. At high temperatures only the layers experi¬ 
encing the external stress move appreciably, but as temperature, T, is 
lowered and the cooperativity becomes pronounced, the perturbation 
is transmitted deeper into the material, reflecting increasing rigidity. 
Note that the extent to which the layers move as depicted, have been 
greatly exaggerated for clarity. 


low the perturbation before fully decaying. Intuitively, one 
would expect that a high temperature liquid, having relatively 
low viscosity, would respond to the external force in a man¬ 
ner such that only the forced layer experiences a substantial 
displacement relative to the opposite boundary. As the tem¬ 
perature is lowered and the viscosity increases, one expects an 
associated increase in the liquids effective, short-lived rigidity. 
This, we argue, corresponds with increasing structural order¬ 
ing and kinetic cohesion of network-like structures in the liq¬ 
uid. At moderate supercooling, then, one expects a deeper 
penetration of shear perturbations and associated displace¬ 
ments sustained by layers of the liquid that are increasingly 
distant from the applied perturbation. As depicted in Figure 
1, the shear is applied to the top layer of the simulation box, 
transverse to the vertical (z) direction. The penetration depth 
is defined as the distance (along the z axis) up to which ap¬ 
preciable effects of shear are observed. At Tg when solid-like 
rigidity has set in, one expects that the whole block of mate¬ 
rial will roughly slide together, such that the penetration depth 
is the length of the material. This is consistent with results 
that show the continuous emergence of a finite shear modu¬ 
lus for temperatures below Tg 1491 . This process is pictured, 
schematically, in Figure Q- Ultimately this length scale is 
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agnostic to the specific type of structural ordering, but can be 
related to cooperativity and the idea of a divergent correlation 
length in ordinary critical theory. 

Models and Methods. Molecular dynamics (MD) simula¬ 
tions were employed using the LAMMPS package tSOll . The 
atoms in the simulation evolved under the infiuence of a semi- 
empirical Finnis-Sinclair type Embedded Atom Model poten¬ 
tial created by Mendelv et al ED with periodic boundary 
conditions. The parameters and coefficients associated with 
the potential were fitted using X-ray diffraction data as well 
as enthalpy of mixing values, and volume measurements in 
the liquid state. This potential has been shown to excellently 
reproduce both the high temperature liquid as well as glassy 
states of NiZr 2 ED 

The simulations were run in the NPT ensemble with 
N=5000 atoms and a target external pressure of P=0. Ther¬ 
mostatting and barostatting were employed using a Nose- 
Hoover thermostat, and barostat respectively, and the velocity- 
verlet algorithm was utilized to integrate the equations of mo¬ 
tion. A 5fs timestep was employed. The initial configurations 
were generated randomly and the atoms were then allowed to 
melt and evolve naturally for 0.25 ns at a temperature of 2200 
K to allow for equilibration. The system was then quenched 
to various target temperatures ranging from 300 K up to 1900 
K using a quench rate ofQ= 10^^ K/s. After the quench, the 
system was allowed to evolve unperturbed for an additional 
0.1 ns. The process was repeated, starting from independent 
initial configurations, for Tg as well as all sampled tempera¬ 
tures below Tg and some representative temperatures above. 

As the glass transition is a kinetic phenomenon without a 
thermodynamic driving force, the glass transition temperature 
Tg is not a constant, and weakly depends on cooling rate, and 
external timescale. Therefore, one has to be careful in identi¬ 
fying its precise location. Typically, various thermodynamic 
parameters show a crossover at the glass transition, associ¬ 
ated with falling out of equilbrium (the system loses its trans¬ 
lational degrees of freedom on the timescale of observation). 
One such property that shows a change in behavior at the glass 
transition is the volume. The temperature dependence of the 
volume shows a ’’kink” at the glass transition temperature, Tg 
(the thermal expansion coefficient, a = has a discontinu¬ 
ity), providing an efficient way to determine Tg. The tempera¬ 
ture dependence of the volume of our system during a quench 
to 300 K is depicted in panel (a) of Figure Q. There appears 
to be a subtle kink in the vicinity of T lOOK. This is in 
good agreement with the results in panel (b) of the same fig¬ 
ure, which was produced by the author of the potential in ED, 
as well as previous numerical work performed under similar 
protocols (521. 

In order to assess that our system is behaving as expected 
before applying shear stresses, we examine the behavior of the 
radial distribution functions at various temperatures. These re¬ 
sults are shown in Figure It is clear that system behaves 
as expected as the glass transition is approached, and in com¬ 
parison with ED, we see that the radial distribution functions 
(RDFs) measured in this work retain the overall shape and 
placement of the peaks. The height of the first peak (and be¬ 
havior of the splitting of the first peak), however, is slightly 




FIG. 2. (Color Online). Panel (a): Specific volume as a function of 
temperature for our simulated system. Panel (b): Specific volume as 
a function of temperature. Reproduced from ED 

different from those in ED and we attribute this to the dis¬ 
crepancy in quench rates. 

We modeled the external shear stress by defining a 4 
angstrom-thick layer at the top of our simulation box, and ap¬ 
plying an external force in the x direction on the atoms in this 
layer. In order to avoid fracturing in such a small system size 
at low temperatures, a force value of only 0.2 eV/A was used. 
A stronger force would be expected to make the effects more 
dramatic, but system size limitations did not allow for higher 
values. The force was left active for 100 timesteps to attempt 
to approximate an impulsive kick at the top of the box. No ex¬ 
ternal forces were applied to the bottom of the box. After an 
observation time of 16,000 timesteps, displacement data 
was extracted. 

Measurement Results. 

To quantify the depth of penetration of the shear stress, we 
plot the displacement of each atom in the shear direction (x- 
direction) versus its position in the transverse height dimen- 
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FIG. 3. (Color Online). Radial Distribution Function at a tempera¬ 
ture of (a): T=300 K < Tg, (b): moderate supercooling with T=1100 
K, and (c): above the melting temperature at T=1500 K. 


sion (z-direction). The displacement represents the net move¬ 
ment in the shear direction from the timestep before the ex¬ 
ternal shear stress was applied, up to the observation time Tq 
(as described in methods section). The height of each particle 
corresponds to the vertical layer it is in at the observation time. 
Due to the periodic boundaries, only atoms in the layers from 
z = I to z = T at the observation time were considered. Figure 
0 shows displacement data for four representative tempera¬ 
tures, i.) deep in the glassy phase (300K), ii.) at Tg (TOOK), 
hi.) in the moderately supercooled regime (1300K), and iv.) 
above Tyneit (1500K). Thermal effects tended to produce large 
motions in the height dimension at temperatures above T^eit^ 
but the effects were not large enough to wash out the effect 
of shear penetration except at very high temperatures (1700- 
1900 K) . Figure ^ serves to quantify the impact of thermal 
noise. In panel (a), the standard deviation of displacement (at 
the observation time) in the direction of applied shear is plot¬ 
ted as a function of the temperature. As expected, the thermal 
noise decreases with decreasing temperatures becoming very 
small as Tg is approached. In panel (b), the average magni¬ 
tude of the particle displacements in the height (z) direction is 
plotted as a function of temperature. It is clear from panel (b) 
that large scale thermal motion may play a role in ’’washing 


out” the shear penetration depth at the highest temperatures 
measured. 

The general response function, of our system to the ex¬ 
ternally imposed shear is a function of distance in the z direc¬ 
tion to the imposed shear, temperature T, and observation time 

To; 

n = n(z,T,Toy (i) 

In this work, we chose a constant observation time, r^, and 
varied the temperature to ascertain the penetration depth along 
the z axis. The penetration depth of the applied shear ulti¬ 
mately has some value depending on the temperature. How¬ 
ever, our ability to extract the exact value depends on the 
observation time chosen. For sufficiently short observation 
times, the effects of the externally applied shear cannot pen¬ 
etrate the system at the lowest temperatures (near and below 
Tg). Therefore, the observation time has to be sufficiently long 
to capture the effect at low temperatures. As shown in Figure 
for very long observation times at high temperatues, the 
effects of the external shear will be nill. Therefore, using an 
observation time which is very long would lead one to con¬ 
clude a much deeper penetration depth at high temperatures 
(see Figure [^. Hence, choosing an appropriate observation 
time is important. For a couple representative temperatures 
we investigated the impact of observation time. In each case 
the duration for which the shear force was applied was con¬ 
stant; for this work we wanted to maintain an approximation 
of an impulsive kick to the system. Oscillatory shears have 
been discussed elsewhere. It was observed that at the lowest 
temperatures (TjTa), the results showed little change when 
observation times were changed by factors of two. For high 
temperatures, the observation time plays a more noticeable 
effect. Afte investigation, we found that the observation time 
employed in this work was sufficient to capture the low tem¬ 
perature effect of the penetration, while not losing the high 
temperature impact except at the hightest temperatures stud¬ 
ied. 

Typically, one would expect the displacement response to 
decay exponentially with depth. As pictured in Figure Q, 
at these system sizes, the shear-induced displacement, while 
not very large, is still quite noticeable. As the displacement 
was not extremely large, we applied a linear regression to the 
data rather than an exponential one. These fits are sufficient 
to quantify the penetration depth. Fits to the data were of the 
form, 

^x(z, To) = m^z-^ So ( 2 ) 

where we denote by S, the displacement in the shear direction 
(x-direction) and z the height of the layer (both measured in 
angstroms). We define the rigidity length scale (the penetra¬ 
tion depth) as 


where m is the slope in Eq. 0- For the temperatures noted 
above, the value of m was averaged over multiple runs, and 
this average was used in Eq. ^ for these temperatures. 
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FIG. 4. (Color Online). Typical Data at four [(a): 300 K, (b): 700 K, (c): 1300 K, (d): 1500 K] representative temperatures both above and 
below Tg 700 K). The red lines correspond to the standard deviation at the two boundaries and center of the material. It is noteworthy that 
they are significantly tighter than the data seems to suggest at this level of zoom. The black lines are the lines of fit from which the slope is 
extracted to define the length scale. Note the dramatically changing behavior as T is lowered. 


The temperature dependence of the length scale (^) is 
shown in Figure 0. A dramatic growth of ^ with decreas¬ 
ing temperature is evident. The first notable penetration of 
the shear, beyond the layer to which the force was applied, 
occurs at a temperature marked Ta{^ 2Tg) (531. Below this 
temperature super-Arrhenius growth of the viscosity may be 
anticipated based upon collective effects (53^5^ . The sud¬ 
den, monotonic increase of ^ at temperatures below pro¬ 
vides direct support to earlier numerical studies, which found 
that metallic liquids begin to develop solid-like features once 
they are cooled below Ta 123 El. These solid-like features 
include the breakdown of the Stokes-Einstein relationship, ex¬ 
ponential stretching of the relaxation functions, and the onset 
of cooperative structural rearrangements during the liquids re¬ 
laxation process m. Indeed, Ta does appear to serve as a 
crossover temperature below which the liquid begins to ex¬ 
hibit a substantial rigid response to external forces, though 
it is a local and transient response. Putting all this together, 
our simulation results allow us to predict that marked growth 
of the shear penetration depth may commence at the same 
temperature as the breakdown of the Stokes-Einstein relation 
in real supercooled liquids. The penetration depth increases 
rapidly as the liquid is supercooled toward the glass transi¬ 
tion temperature, Tg. Below the latter temperature, the ma¬ 


terial is glassy and exhibits structural rigidity on all practical 
timescales, such that shear perturbations propagate the length 
of the material and appear to diverge. 

As discussed previously, we performed multiple (typically 
six) independent measurements at all temperatures T < Tg ^ 
lOOK, as well as most of the representative temperatures 
above Tg (300-800, 1300, 1400, 1700K). At each of these 
temperatures, averaging was done to determine the slope, m, 
in the fit of Eq. 0- The apparently periodic nature of ^ be¬ 
low Tg can be attributed to noise due to the lengthscale being 
essentially divergent to the system size at these temperatures. 
This point is vividly made in Eigure ([^, which depicts the 
average slopes, m, as a function of temperature. The error 
bars on the points with multiple runs corresponds to the stan¬ 
dard deviation in slopes. We see that, at temperatures near 
and below Tg, the combination of the average value and as¬ 
sociated error bars lead to m being virtually indistinguishable 
from zero, consistent with the divergence of the penetration 
depth beyond the system size at the glass transition. It should 
also be immediately noticeable that the data points associated 
with the value of the penetration depth at temperatures T=600, 
650, and 750 K, are not plotted in Eigure 0. This is because, 
as seen in Eigure the signal to noise ratio (^) for these 
points was a factor of three for the data point at 600K and a 





































6 


lU" 

Q _ 

(a) 

♦ 

o 

Oxih ^ 

A. - 


♦ 

4 

T 

1 

G 

1 .. 

♦ 


2 

0 

: 


♦ 


; ^ 

. 

u* 

^- 

r . 


0 500 1000 1500 2000 


8 
6 

I^(A) 4 
2 

0 - 

0 500 1000 1500 2000 

T(K) 

FIG. 5. (Color Online). Panel (a): Plot of standard deviation of 
displacement in direction (x-axis) of applied shear as a function of 
temperature. Panel (b): Average displacement in height dimension 
(z-dimension) as a function of the temperature. 
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FIG. 6. (Color Online) Displacement data for a random configu¬ 
ration at 1700 K. Panel (b) is the displacement for the observation 
time (16,000 timesteps), panel (a) is the data at an earlier time (3000 
timesteps after the shear is turned off). As panel (b) shows, at high 
temperature the shear induced displacements appear to be far smaller 
at the standard observation time used in this work. Nevertheless, at 
earlier times, as seen in panel (a), the displacements are much more 
noticeable. 


factor of ten for the data points at 650 and 750 K. This large 
relative error is due to the fact that the length scale, is ex¬ 
ceedingly large as the corresponding average slope m is very 
small, and in fact virtually indistinguishable from zero (see 
Eq. 0). This corresponds to total penetration of the shear to 
beyond the system size, and the fluctuations about m=0, are 
to be expected due to ordinary thermal effects. Because of 
the large relative error in the aformentioned data points, we 
removed these data points from Figure 0 so as to not mask 
the overall monotonic increase of ^ with incredible values. It 
may at first seem concerning that the data point at T=750 K is 
suggestive of near divergene considering it is above Tg. This 
is, in fact, not an issue, as the value of Tg is not precise, and 
very likely falls within the T=700 to 750 K range for this sys¬ 
tem size and quench rate. Also, the close proximitiy to Tg and 
limitations of resolution at this system size, would lead to the 
impact of the arrest at the glass transition strongly influenc¬ 
ing temperatures asymptotically close to Tg. Clearly, Figure 
^ serves not only to explain the fluctuations but also rein¬ 
forces the idea that the penetration depth appears to diverge to 
the system size in the vicinity of Tg, and is perhaps the most 
consequential and rigorous result in this work. 

As a lower bound on the shear penetration depth, in Figure 
0 we plot 


Because this is a lower bound, we can strongly assert, based 


on our data, that the penetration does indeed show dramatic 
increase upon supercooling. For temperatures below Tg, the 
length scale is so large and slopes so small, that the observed 
fluctuations of ^ may be statistical (see Figure [^. We also 
conclude that the length scale becomes, at least, considerably 
larger than the system size at Tg, and may in fact diverge. For 
the impulsive kick we applied this will remain the case at all 
longer timescales of observation. However, if one were to ap¬ 
ply a static shear at the top of the box and left this shear on for 
a time longer than the relaxation time, then even below Tg, the 
length may not diverge. On all practical timescales, though, 
it would, not diminishing this length as a natural candidate 
for the glass transition problem. The precise behavior of the 
penetration depth as the duration of the static shear stress is 
varied is an interesting problem, but requires a different type 
of analysis, and will be addressed in a future work. 

Scaling Arguments. Previous studies have examined the be¬ 
havior of various proposed length scales in the vicinity of Tg 
(ortheVogel-Fulcher-Tammann temperature To ir7l l57ti59l . In 
(301, a diverging length scale associated with liquid-like de¬ 
fects at Tg produced an exponent y=l; This value constitutes 
an upper bound on the exponents reported in other works. 
Researchers in (601 found a critical exponent of y=0.875 
for the largest icosahedral cluster size in a model metallic 
glass former (56l. In both ED, focusing on inherent struc¬ 
tures in a binary Lennard-Jones glass former, and in (62l, 
by largely studying medium-range bond orientational order in 
colloidal liquids, scaling analyses gave exponents of y 2/3, 
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FIG. 7. (Color Online). Plot of the length scale, ^ versus temperature. 
All temperatures below Tg were averaged over multiple independent 
runs, as were select, representative temperatures above Tg. 
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FIG. 8. (Color Online). Plot of the average value of the slopes, 
m, at each of the measured temperatures along with their associated 
standard deviations cr (when multiple runs were performed). 
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FIG. 9. (Color Online). Lower bound, ^lb, on the shear penetration 
depth. See text. 


FIG. 10. (Color Online). Power law fits to the shear penetration 
depth as a function of reduced temperature (measured relative to 

glass transition temperature Tg). Panel (a): ^ oc • Panel 



in rough agreement with a three-dimensional Ising exponent 
(y=0.625(1) 16^ ). We performed a similar scaling analysis of 
our data, fitting a function of the form 


^oc 



(5) 


To represent the value of the length scale at Tg itself, we 
interpolated the value at T=710K using the line connecting 
T=800K and T=700K. In addition, the data point at T=750 K, 
was excluded from this scaling analysis for reasons discussed 
above. When using the full range of temperatures in applying 
the power law fit, we extracted an exponent of y=l. This is 
depicted in Figure ( p^ , panel (a). Decreasing the temperature 
range considered in the scaling to only include temperatures 
very close to Tg, caused the value of the exponent y to de¬ 
crease. This is depicted in panel (b) of Figure ( p^ , where a 
value of y=0.713 was found. Clearly, as the scaling is applied 
to a more and more asymptotic region around Tg, the value of 
y appears to approach a value consistent with Ising-like scal¬ 
ing. Our value for , is, thus, in rough agreement with pre¬ 
viously suggested exponents extracted by different means. A 
similar study conducted in d studied the high temperature 
(T > Tm) correlations of the anisotropic part of the atomic 
level stress. In this work the authors found an exponent of 
y ^Q.l for the low temperature extrapolation of the correla- 
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FIG. 11. (Color Online). Power law scaling, ^ oc in the 

asymptotic region below the crossover temperature, ra¬ 


tions in their two-dimensional system (the high crystallization 
rate at temperatures below thwarted a direct study at low 
temperatures). Our possible scaling may provide further evi¬ 
dence of a universal nature of the length scale at deep super¬ 
cooling. As Tg is not a true thermodynamic temperature, it is 
unclear what this scaling may mean, but it may be suggestive 
of universality in the glass transition. 

In [69, 70] theoretical arguments for the scaling of the 
length scale about the crossover temperature, Ta, were pro¬ 
vided. It was suggested when asymptotically approaching 
Ta from below, that a characteristic structural domain size, 
1, scaled as 


1 = ta^\ 


( 6 ) 


where ta = . We examined the shear penetration depth 

as a function of the reduced temperature, ta, as depicted in 


Figure ( pT] ). We observed that a power law is a good fit in this 
region with an exponent, va ~ 1.335. If the lowest temperature 
(highest reduced temperature) point in Figure ( pTj ) is removed, 
then a value of va ~1.5 will be obtained instead. 

Conclusion. We have defined the shear penetration depth as 
the distance over which supercooled liquids can support shear 
appreciably. This definition is based on a simple physical pic¬ 
ture of liquids as a continuum unable to globally support shear 
stress. With decreasing temperature the viscous inter-layer 
forces, rigidity, and the lifetime of connectivity increase. At 
Tg, the glass transition temperature, these quantities mirror 
those of crystals and the material is solid. This scale has the 
added benefit that it may be readily experimentally accessi¬ 
ble. Results from simulations on NiZr 2 , a typical fragile glass, 
support this definition. In the current work, we applied, in 
simulatum, a shear force to the top layer of a NiZr 2 system. 



FIG. 12. (Color Online). Depiction of growing interconnectivity of 
icosahedral clusters with supercooling in a very similar 

metallic glass former. The Cu atoms are marked red and the Zr by 
purple. In Panel (a) we show the longest interconnected cluster at 
1200 K. Panel (b) shows the longest connected cluster at Tg (800 K 
for this system). Note that interconnecting icosahedra percolate at 
T^. (These results are similar to those in ED.) 


By measuring the penetration depth, we have demonstrated 
a dramatically increasing structural scale upon supercooling 
toward Tg. This leads to the conclusion that the shear penetra¬ 
tion depth marks a very natural candidate for the structural 
length scale characterizing glassy dynamics. Furthermore, 
and of equal importance, the shear penetration depth can be 
measured experimentally. While it might be practically diffi¬ 
cult, in theory, the penetration depth can accessed experimen¬ 
tally in a way that does not rely solely on scattering, and may, 
therefore, be easier to investigate. This is a major advanatage 
for this lengthscale, and sets it apart from previously proposed 
lengths. 

We ultimately believe that the shear penetration depth is 
intimately connected to the structure of the supercooled liq¬ 
uid. It has been shown 15411561164^691 in extensive numerical 
studies that clusters of locally preferred structural order tend 
to grow and interconnect as temperature is lowered in super¬ 
cooled liquids. These clusters locally minimize the relevant 
free energy and hence are stronger and more stable to fluctu¬ 
ations. Interconnections of the clusters increase in length and 
lifetimes HaisslIMl upon lowering temperature, and even¬ 
tually span the system size at the glass transition. This per- 
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eolation of a structural network, which is depicted in Figure 
( p^ for a cousin configuration, is very likely the source of the 
penetration depth as well as a leading cause of the arrest at the 
glass transition temperature Tg. 

Due to their locally stable nature and tight binding, the 
clusters resist thermal breakup and lock into a rigid structure 
forming a force network that can propagate shear. The inter¬ 
locking and cohesiveness of this network also serves to slow 
down the dynamics (661 . as sufficiently large thermal fluctu¬ 
ations are needed to break the network and this becomes less 
likely with lower temperature. In fragile glasses, this network 
has to form quickly over the temperature range encountered 
in typical experiments. This is likely due to the largely non- 
directional binding in fragile glasses which lacks the natural 
network found in strong covalent liquids. 

In metallic liquids, the network is likely icosahedral (541 [5^ 
[64l[65][67l. In silicates (typical of the strong classification) a 
natural tetrahedral network with strong bonds and direction¬ 
ality is present. It has been suggested that (67l networks of 
locally preferred structures tend to form in fragile glasses be¬ 
ing either icosahedral or crystal-like at short range. 

The notion that the shear penetration in amorphous solids is 
due to a system spanning network is a universal one. As dis¬ 
cussed a network forms over a narrow range in fragile glasses 
leading to the super-Arrhenius increase of viscosity and caus¬ 


ing the rigidity. In strong glasses, a tetrahedral network forms 
at high temperature and becomes increasingly cohesive as the 
temperature is lowered to Tg. Other forms of a system span¬ 
ning network can also exist. In colloids a frictional or con¬ 
tact network can be created by jamming, and in fact a rigidity 
length scale has been proposed for these systems (26l. The 
formation of a contact network has also been shown to occur, 
albeit short-lived, in some discontinuous shear thickening flu¬ 
ids (68]| . This contact network may also play a role strong to 
fragile crossovers in high pressure thermal glasses. 

Based on the above discussion, it is clear that a shear pene¬ 
tration length scale can be quite naturally extended to many, if 
not most, glassy systems. This leads naturally to the connec¬ 
tion between slowing down and network formation. The fact 
that our length scale begins to grow substantially, only when 
supercooled beneath Ta, further suggests a structural origin 
for the glass transition. 
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